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Abstract 

Pesticides have the potential to leave harmful effects on humans, animals, other living organisms, and the en- 
vironment. Several human metabolic proteins inhibited after exposure to organophosphorus pesticides absorbed 
through the skin, inhalation, eyes and oral mucosa, are most important targets for this interaction study. The crys- 
tal structure of five different proteins, PDBIDs: 3LII, 3NXU, 4GTU, 2XJ1 and 1YXA in Homo sapiens (H. sa- 
piens), interact with organophosphorus pesticides at the molecular level. The 3-D structures were found to be of 
good quality and validated through PROCHECK, ERRAT and ProSA servers. The results show that the binding 
energy is maximum -45.21 relative units of cytochrome P450 protein with phosmet pesticide. In terms of H-bond- 
ing, methyl parathion and parathion with acetylcholinesterase protein, parathion, methylparathion and phosmet 
with protein kinase C show the highest interaction. We conclude that these organophosphorus pesticides are more 
toxic and inhibit enzymatic activity by interrupting the metabolic pathways in H. sapiens. 

Keywords: docking, organophosphorus pesticides, comparative analysis, modeling, toxicity analysis, acetylcho- 
linesterases, Homo sapiens 



INTRODUCTION 

Pesticides, as a consequence of massive use in ag- 
riculture and other human activities, are widely dis- 
tributed environmental contaminants and are subject 
to restricted use by legislations aimed at the protec- 
tion of natural ecosystems and human health in India, 
Europe, the USA, and many other countries in the 
world. A wide structural variability also characterizes 
the pesticide subfamilies like insecticides, herbicides 
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and fungicides, that are grouped together according to 
the target of biocide activity [1] . Pesticides, herbicides 
and insecticides have the potential to exert harmful 
effects on humans, animals, birds and other living or- 
ganisms via inhalation or breathing, skin contact and 
food consumption. They also destroy our biodiversity 
and environment. Some possible changes occur after 
exposure such as skin irritations, nausea, breathing 
problems, cancer, hematological disorders, damage to 
the reproductive organs, and neurological disorders. 
The different routes of entry and the high binding ef- 
ficiency of presticides with particular targets proteins 
in the human body have been taken into account. 

Organophosphorus pesticides are the most fre- 
quently used pesticides in the world, with uses ranging 
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from commercial and home use to agricultural appli- 
cations for controlling unwanted insect pests [2] . Orga- 
nophosphate pesticides include parathion, malathion, 
methylparathion, chlorpyrifos, diazinon, dichlorvos, 
phosmet, monocrotophos, fenthion, quinalphos, tetra- 
chlorvinphos and azinphosmethyl. Organophosphorus 
pesticides obtain their toxicity from their ability to in- 
hibit many metabolic and physiological enzymes like 
acetylcholinesterase (AchE), cytochrome P450, pro- 
tein kinase C, and glutathione S -transferases (GSTs), 
causing neurotoxicity in humans [3] . For example, the 
presence of AchE in insects, birds, fish and all mam- 
mals give this class of pesticides enormous toxicity 
towards unintended targets to disrupt the endocrine, 
metabolic and digestive systems in the human body [3] . 

Human AchEs are a class of enzymes which cata- 
lyze the hydrolysis of acetylcholine (Ach), an ester 
which acts as a neurotransmitter. The reaction cata- 
lyzed by AchE is: Ach + H 2 0 choline + acetate [3] . 
The inhibition of AchE by organophosphorus pesti- 
cides occurs as a result of the phosphorylation of the 
serine residue in the active site. The hydroxyl group 
of serine residue acts as an electrophile which attacks 
the nucleophilic phosphorus. After phosphorylation, 
the protein is highly stable and the hydrolysis of Ach 
is blocked. In some cases, depending on the chemi- 
cal structure of the pesticide, the phosphorylation and 
inhibition may be irreversible [3] . The inhibition of the 
protein causes formation of Ach in the neural synaps- 
es. The main function of AchE is to recycle Ach by its 
hydrolysis at cholinergic synapses in order to restore 
the membrane potential after propagation of a nerve 
impulse [3] . The principle of cholinesterase inhibition 
i.e. in a biosensor is used as a means of detection of a 
phosphorothionate ester compound and inhibition of 
AchE may be low or absent [4] . It is reported that AchE 
is a target enzyme for biologically active compounds 
ranging from anti-Alzheimer disease agents acting as 
reversible inhibitors to organophosphorus pesticides [5] 
and warfare agents which act as reversible or irrevers- 
ible inhibitors. 

The activation of organophosphate pesticides has 
been attributed to the cytochrome P450 family of en- 
zymes. P450s have also been shown to carry out direct 
detoxification of organophosphorus pesticides through 
dearylation. Humans have the ability to activate or- 
ganophosphorus pesticides through CYP1A2, 2B6, 
2C19, and 3A4 cytochromes, which are particularly 
sensitive to their actions [6] . This hypothesis reported 
that the activation and detoxification of parathion 
and chlorpyrifos pesticides are found in human liver 
microsomes [7] . In addition, the recombinant human 
P450 protein was used to quantify organophosphorus 



pesticides by human liver microsomes and the kinetic 
values may vary widely because of the marked vari- 
ability of P450 content and activity in procured human 
specimens and differences in incubation conditions 
and analytical methods [8] . 

The third protein which we target are GSTs, a 
widely distributed family of detoxifying dimeric 
enzymes found in most forms of life. GSTs inacti- 
vate toxic effects by chemically bonding them to the 
tripeptide glutathione, making them soluble so that 
the body can easily excrete them. Pathogenic parasites 
also make their own GSTs to help them inactivate the 
drugs [9] . Resistance to organophosphorus pesticides is 
considered to be due to the metabolism of these com- 
pounds by GSTs [10] . Many researchers state that in- 
secticide resistant insects have elevated levels of GST 
activity in crude homogenates, which suggests a role 
for GSTs in resistance [10] . Multiple forms of these en- 
zymes have been reported for mosquitoes, house fly, 
Drosophila, sheep blow fly and grass grub [nl . GSTs 
are also involved in intracellular transport, biosyn- 
thesis of hormones and protection against oxidative 
stress. In addition, they contribute to the removal of 
toxic oxygen free radical species produced through 
the action of pesticides. They have peroxidase [12] and 
isomerase activity [13] , which inhibits the junction of 
N-terminal kinase (thus protecting cells against H 2 0 2 - 
induced cell death) and they are able to non-catalyti- 
cally bind a wide range of endogenous and exogenous 
ligands [14] . 

The fourth protein, i.e. protein kinase C (PKC), is 
also one of the most key target proteins for organo- 
phosphate pesticide toxicity, because it can also have 
an effect on signaling pathways by activating PKC [15] . 
PKC may mediate signaling effects and to date there 
has been no information regarding whether these xe- 
nobiotics can modulate PKC, which is a significant 
event signaling the increase in endothelial perme- 
ability and cell proliferation. However, the activation 
is probably indirect through organophosphorus pesti- 
cides-mediated formation of reactive oxygen species 
(ROS) that activate PKC, which can be measured in 
human liver and brain cytosolic fractions [15] . Organo- 
phosphorus pesticides will induce ROS formation and 
oxidative stress has also been shown to be associated 
with apoptosis in different tissues. Organophosphorus 
pesticides can also inhibit the steroid androgen recep- 
tor (AR), which causes steroid hormone disturbances 
in the human body. It is reported that in hepatic tis- 
sues, the greatest increase in activities was observed 
with TCDD (2, 3, 7, 8-tetrachlorodibenzo-p-dioxin; 
herbicide), chlorpyrifos, endrin and Cd (II), while 
chlorpyrifos and fenthion exerted the greatest increas- 
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es in the brain tissues . Protein kinase C may be a 
best target protein of free radicals and oxidative stress, 
leading to altered cell proliferation and differentia- 

[15] 

tion . 

The last protein is al-antichymotrypsin (ACT), an 
inhibitor of proteinases of the chymotrypsin class. It is 
also activated in Alzheimer's disease patients. ACT is 
synthesized primarily in hepatocytes and secreted into 
the blood [16] . The expression of ACT in hepatic cells is 
known to be enhanced by interleukin-6 (IL-6), to some 
extent by IL-1, and also by glucocorticoids^ 73 . ACT is 
synthesized in human bronchial and breast epithelial 
cells, epididymal cells, predominantly in the choroid 
plexus in normal brain, and in astrocytes or astroglia, 
and to a small extent in monocytes [16] . Brain micro- 
vessel endothelial cells also release ACT [17] , indicating 
that the inhibitor may also perform unique functions in 
local microenvironments. ACT may be also involved 
in controlling the oxidative damage because of its 
correlation to inhibition of oxygen consumption and 
superoxide generation in human granulocytes. Com- 
plexes formed by ACT and chymotrypsin regulate the 
production of superoxides in neutrophil membranes 
by interacting with NADHP oxidase [18] . 

Molecular docking is a frequently used method in 
computer-aided drug design. It evaluates how small 
molecules called ligands like organophosphate pes- 
ticides and the target macromolecules (e.g. receptor, 
enzyme or nucleic acid) fit together. Hence, this re- 
search hypothesis focuses on a comparison and in- 
teraction study of five different metabolic and physi- 
ological proteins [PDBIDs: 3LII, 3NXU, 4GTU, 2XJ1 
and 1YXA in Homo sapiens (H. sapiens)], which 
are inhibited by the exposure to ten different organo- 
phosphorus pesticides and exert their toxic effects on 
humans metabolism, by using online bioinformatics 
tools and softwares. We are paying attention to the 
interaction of organophosphorus pesticides with target 
proteins and predict which one is more toxic and re- 
veal the injurious effects on human by dermal layers 
exposure, inhalation, eyes exposure and oral contact. 

MATERIALS AND METHODS 

Model description 

The crystal structures of five different proteins, 
AchE, cytochrome P450, GST, PKC and ACT, were 
obtained from the Protein Data Bank (PDB) website 
(http://www.rcsb.org/) in *.pdb (dot pdb) format [19] . 
This format is recommended for protein structures 
that are obtained by X-ray crystallography or NMR 
studies. For further analysis, we selected recombinant 
human AchE (PDB ID: 3LII) [20] , which has a 3.20 A 



resolution. The second one was the crystal structure 
of human cytochrome P450 3 A4 bound to an inhibitor 
ritonavir (PDB ID: 3NXU) [21] and it has a 2.0 A reso- 
lution. The third was ligand-free homo dimeric human 
GST m4-4 (PDB ID: 4GTU) [22] , which has a 3.30 A 
resolution. The protein kinase pim-1 in complex with 
a small molecule inhibitor (PDB ID: 2XJ1) [23] has a 
2.13 A resolution and serpina3n, a murine orthologue 
of human al-antichymotrypsin (PDB ID: 1YXA) [24] 
has a 2.10 A resolution. Crystallographic waters and 
hetero atoms were removed and the structures were 
fully solvated before docking analysis. 

Sequence alignment 

AchE, P450, GST, PKC and ACT are composed of 
534, 457, 217, 273, and 372 amino acids, respectively. 
The protein sequence files were obtained from the 
PDB in *.txt (dot txt) format. Multiple alignments of 
the related sequences were performed using the Clus- 
talW program accessible through the European Bio- 
informatics Institute website (http://www.ebi.ac.uk/ 
Tools/clustalw2/index.html) [25] . 

Secondary structure prediction 

Secondary structure analyses of AchE, P450, GST, 
PKC and ACT in humans were described in Protein 
Data Bank. Prediction of different domains of all these 
proteins was carried out through the SUPERFAMILY 
sequence search server [26] by using the hidden Markov 
models. 

Model validation 

The validation of AchE, P450, GST, PKC and ACT 
protein structures were assessed by using Ramachan- 
dran plot through PROCHECK validation package [27] . 
The initial model showed problems in the conforma- 
tion of different loop regions, which were subjected 
for loop modeling and further validated by ERRAT 
web server [28] . ERRAT plot gives measurement of the 
structural error for each residue in the 3-D structure 
model. This process was repeated iteratively until 
most of the amino acid residues were below 95% cut- 
off value and the residues which lie above 95% cut- 
off value were subjected to loop modeling in MOD- 
ELLER. Finally, all protein models showing the best 
PROCHECK and ERRAT plot were subjected to na- 
tive protein folding energy evaluation by using ProSA 
program [29] . ProSA is an interactive web which re- 
quires the atomic coordinates of the model to be eval- 
uated and recognition of errors in three-dimensional 
structure. 

Protein function analysis 

We also analyzed the novel functions of AchE, 
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P450, GST, PKC and ACT proteins by means of 
SVMProt server (BIDD server) with the aim of Sup- 
port Vector Machine (SVM) learning techniques, 
which classify proteins into functional families from 

[30] 

its primary sequences . 

Ligand binding site (active site) prediction 

Pocket-Finder or Q-site finder is a molecule-bind- 
ing site prediction server based on the Ligsite algo- 
rithm and also compares with the CASTp server 131 ' 32] . 
It works by scanning a probe radius 1.6 A along all 
gridlines at a grid resolution of 0.9 A surrounding the 
protein. The probe also scans cubic diagonals. Grid 
points are defined to be part of a site when the probe 
is within range of protein atoms followed by free 
space followed by protein atoms. Grid points are only 
retained if they are defined to be part of a site at least 
five times. 

Protein-ligand interaction study (docking) 

The 3-D structures of AchE, P450, GST, PKC and 
ACT were further used for insilco docking study to 
know the interaction between the organophosphorus 
pesticides and target proteins. Various methods ap- 
plied in this study are given below 

Preparation of proteins 

The crystal structures of AchE, P450, GST, PKC, 
and ACT and their PDBID, 3LII, 3NXU, 4GTU, 
2XJ1, and 1YXA [20 24] were downloaded from the 
PDB in *.pdb (dot pdb) format. The crystal structures 
of proteins have already presented some ligand mol- 
ecules, which were removed. The other necessary ac- 
tion was to remove water molecules from the surface 
of the protein molecules. It is compulsory because the 
extra water molecules will mask the protein surface 
from the ligand. All protein 3-D models were prepared 
for docking by removing waters and hetero atoms for 
docking analysis by using PyMol V2.0.7 software [33] . 

Preparation of ligand 

The three dimensional structures of ten orga- 
nophosphate pesticides were downloaded from the 
Pubchem compound search of NCBI database [34] . The 
structures were downloaded as *.sdf (dot sdf) file 
format. The 3-D structure of ten organophosphate 
pesticides like malathion, parathion, methylparathion, 
monocrotophos, chlorpyrifos, fenthion, quinalphos, 
phosmet, tetrachlorvinphos and azinphosmethyl have 
Pubchem IDs are: CID4004, CID991, CID4130, 
CID5371562, CID2730, CID3346, CID26124, 
CID12901, CID5284462 and CID2268, respectively. 
The 3-D structures of organophosphate pesticides 
were open through PyMol V2.0.7 [33] , a visualization 



tool and converted into *.pdb file. 

Docking studies using PatchDock and FireDock 

The crystal structures of AchE, P450, GST, PKC 
and ACT were used for docking analysis through 
PatchDock (http://bioinfo3d.cs.tau.ac.il/), where can- 
didate solutions were generated by rigid-body docking 
methods [35] . PatchDock determined the best starting 
candidate solutions based on shape complementarily 
of soft molecular surfaces of proteins. The Clustering 
RMSD was 4.0 A for analysis and complex type was 
set to as default. The PatchDock algorithm divides the 
Connolly dot surface representation of the molecules 
into concave, convex and flat patches. Then, com- 
plementary patches are matched in order to generate 
candidate transformations [35] . Each candidate trans- 
formation is further evaluated by scoring function that 
considers both geometric fit and atomic desolvation 
energy. The 1000 best docked candidate transforms 
from PatchDock, based on global energy, attractive 
and repulsive van der Wall's interactions, partial elec- 
trostatics, atomic contact energy (ACE), and addition- 
al estimations of the binding free energy were used in 
FireDock (http://bioinfo3d.cs.tau.ac.il/) [36] . FireDock 
re-scored the 10 top candidate solutions by restricting 
the flexibility to the side-chains of the interacting sur- 
face and allowing small rigid-body movements. This 
study was carried out by selecting the first best can- 
didate solution from FireDock were retained and then 
visualized through PyMol V2.0.7 [33] . 

RESULTS 

Model description 

The length of amino acid residues of AchE, P450, 
GST, PKC and ACT varies between 217 and 534 
amino acids. The downloaded structures are described 
as text files, which contain necessary information on 
the molecule such as the number of atoms, name of 
atoms, bond distances, angles, dihedral angles, and the 
number of residues. Multiple sequence alignments of 
five amino acid sequences were carried out through 
the ClustalW program, which showed that sequences 
were not identical to each other. 

The three dimensional model of five proteins in H. 
sapiens were validated by VERIFY 3-D score pre- 
dicted by the SAVES server [37] . VERIFY 3-D analyzes 
the compatibility of an atomic model (3-D) with its 
own amino acid sequence. For all 3-D models, the 
range varies between 94.10% and 98.54% for ACT 
and PKC, respectively. The residues have a score of 
greater than 0.2, which indicates a good quality model 
(Table 1). The ERRAT server is used for analyzing the 
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statistics of non-bonded interactions between differ- 
ent atom types and scores greater than 50 are normally 
acceptable. For all five protein models, ERRAT score 
varies between 98.256 and 87.833 of 1YXA (ACT) 
and 3LII (AchE), respectively, which fall within nor- 
mal range for high quality models (Fig. 1 and Table 1). 
ProSA is widely used to check 3-D models of protein 
structures for potential errors. The Z-score indicates 
overall model quality and measures the deviation of 
the total energy of structure with respect to an en- 



ergy distribution derived from random conformations. 
The ProSA score was negative for the modeled pro- 
tein, which indicates its validity. The ProSA profiles 
calculated the protein structures, which were found 
similar to the energy of all five different structures of 
protein from PDB listed in Table 1. The overall model 
quality of all five 3-D structures was reflected by the 
minimum Z-score, which was -7.63 for PKC (PDBID: 
2XJ1) and the maximum Z-score, which was -10.62 
for AchE (PDBID: 3LII) as is shown in Table 1. 



Table 1 Ramachandran map of AChE, P450, GST, PKC and ACT proteins there PDBIDs are 3LII, 3NXU, 4GTU, 
2XJ1 and 1YXA respectively were calculated with the PROCHECK program 



Model Description 


3LII 


3NXU 


4GTU 


2XJ1 


1YXA 


Residues in most favoured regions [A,B,L] 


358/82.1% 


356/90.4% 


217/91.4% 


224/94.1% 


308/91.7% 


Residues in additional allowed regions [a,b,l,p] 


71/16.3% 


36/9.1% 


8/4.0% 


12/5.0% 


28/8.3% 


Residues in generously allowed regions [~a,~b,~l,~p] 


4/0.9% 


1/0.3% 


1/0.6% 


1/0.4% 


0/0.0% 


Residues in disallowed regions 


3/0.7% 


1/0.3% 


3/0.7% 


1/0.4% 


0/0.0% 


ERRAT overall quality score 


87.833 


97.105 


90.431 


97.619 


98.256 


VERIFY 3-D (3D-1D score > 0.2) 


97.76% 


96.29% 


95.54% 


98.54% 


94.10% 


ProSA Z-score 


-10.62 


-9.84 


-8.93 


-7.63 


-8.64 


SCOP (Domains) 


Region: 3-530 


Region: 5-457 


Region: 85-216 


Region: 4-263 


Region: 1-371 








Region: 2-84 







The Z-score and ERRAT overall quality score by servers like ERRAT and SUPERFAMILY sequence search server. 



Model validation 

The geometry of AchE, P450, GST, PKC and ACT 
3-D structures were evaluated through Ramachan- 
dran plot calculations by using PROCHECK. Stereo- 
chemical evaluation of backbone Psi and Phi (O) 
dihedral angles of five human proteins were revealed 
in different percentages i.e. 82.1%-94.1%, 4.0%- 
16.3% and l%-4% residues were diminishing within 
the most favored regions, additionally allowed regions 
and generously allowed regions, respectively (Table 
1). The dihedral angles revealed that some residues 
like l%-3% disallowed regions of Ramachandran plot. 
The model has normal distribution of residue types 
over the inside and the outside of the protein struc- 
tures. The residues in the disallowed region were ig- 
nored as they were not present near the active site nor 
were they involved in ligand binding. 

By using the SUPERFAMILY sequence search 
server, the crystal structure of AchE consists of one 
interacting domain that belongs to the family of 
AchE, an alpha/beta-hydrolase N-terminal domain 
(residues 5-530) and their E-value is 4.34e-163. The 
cytochrome P450 (3NXU) has one domain that lies 
between the 5 and 457 amino acids and belongs to 
the cytochrome P450 superfamily and the E-value is 
1.83e-135. The third one is GST (4GTU), which has 
two domains that lie between region 85-216 and 2-84 
and their superfamily are GST C-terminal domain-like 



and thioredoxin and the expected E-value was 1.77e-46 
and 1.67e-23, respectively. Similarly, in the fourth and 
fifth proteins, PKC (2XJ1) and ACT (1YXA) have 
one family that lies between amino acids 4-263 and 
1-371 and belongs to protein kinase-like (PK-like) and 
serpins superfamily, and their E-value was 1.25e-71 
and 2.23e-135, respectively, that is assigned by SCOP 
domains. 

Prediction of secondary structures 

The secondary structure prediction, which was 
updated in PDB used in the analysis, revealed that 
random coils dominated among secondary structure 
elements followed by alpha helix, extended strand and 
beta turns. The crystal structure of AchE, P450, GST, 
PKC and ACT consists of a total of 8 to 21 a-helices 
in GST and P450 proteins and the number of p-sheets 
varies between 4-16 in GST and AchE, respectively, 
as is shown in Fig. 2. The N-terminal domain consists 
of two anti-parallel p-sheets, forwarding two P barrel 
domains in 3LII and 3NXU crystal structure. The hu- 
man AchE protein (3LII), which consists of 16 mixed 
p-sheets and 20 large and small a-helices, similarly 
in P450 (3NXU), has one anti-parallel p-sheet, seven 
mixed p-sheets and 21 a-helices. The third one is GST 
(4GTU), which has 2 pairs of anti-parallel p-sheets 
and eight a-helices, and the fourth is the PKC (2XJ1), 
which has 12 mixed p-sheets and 13 a-helices. The 
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Fig. 1 The overall quality score (E-value) of A, B, C, D, E represents AchE, P450, GST, PKC and ACT proteins 
respectively by using ERRAT server. 



last is ACT (1YXA), which contains 14 p-sandwiches 
and 12 a-helices (Fig. 2). 

Functional analysis by support vector machine 

Different unknown and hidden functions of five 
different human proteins were predicted using ma- 
chine learning technique like a statistical support vec- 
tor machine-based classifier i.e. SVMProt (Fig. 3). 
The comparative analysis of AchE, P450, GST, PKC 
and ACT for functional assignment shows that it be- 
longs to the transferase group of proteins as shown in 
Fig. 3. AchE (3LII) and P450 (3NXU) belong to the 
transmembrane region proteins. Human proteins like 
AchE, PKC and al-antichymotrypsin have magnesium 
binding (58.6%), metal-binding properties. AchE and 
P450 have the tendency to bind with iron (97.0% and 



99.1%) and also copper (58.6%) properties, respec- 
tively. Comparative analyses have shown that differ- 
ent domains of al-antichymotrypsin belong to ATP- 
binding cassette (ABC) function shown in Fig. 3. 

Ligand binding site analysis 

The potential ligand binding sites (LBSs) of AchE, 
P450, GST, PKC and ACT were identified by Pocket 
Finder program. A total of ten possible binding sites 
were obtained in all five proteins. The possible ac- 
tive sites obtained from the CASTp server in all 
five proteins are shown in Fig. 4. The frequently 
involved amino acid residues in human AchE are 
Gln71, Tyr72, Asp74, Glu81, Thr 83, Tyrl24, 
Glu202, Ser203, Trp236, Tyr337, Pro368, Arg463 
and Asn533. The involved amino acid residues in 
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Fig. 2 Secondary structure of five different human proteins. A: human acetylcholinesterase (3LII). B: cytochromeP450 
(3NXU). C: gutathione S-transferase (4GTU). D: represents protein kinase C (2XJ1). E: is alpha- 1-antichymotryp sin (1YXA). The 
arrow indicates the P-sheets with flexible residues shown in yellow, a-helices in red color and green are the loop regions. The protein 
structures were visualized through PyMol V2.0.7. 
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Protein Family Name 


EL 

Value 


Value 

TO 


Iron-binding 


3.7 


97.0 


Transmembrane 


3.5 


96.4 


Aptamer-binding protein 


3.2 


95.2 


Antigen 


1.5 


73.2 


Structural protein (Matrix protein, 
Core protein, Viral occlusion 
body, Keratin) 


1.5 
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Metal-binding 


1.3 


6%. 5 


EC 4.2.-.-: Lyases - Carbon- 
Oxygen Lyases 


1.3 
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Photosystem I 


1.0 


52.6 


Copper-binding 
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52.6 


Magnesium-binding 
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52.6 



Protein Family 


Value 


P- Value 

(%) 


EC 1.14.-.-: Qxidoreductases 






- Acting on paired donors 
with incorporation or 
reduction of molecular 


7.1 


99 .0 


oxygen 






Transmembrane 
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99 .0 


Iron-binding 
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6.4 


99.1 
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Metal-binding 


2.2 
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EC 2.7.-.-: Transferases - 






Transferring Phosphorus- 
Containing Groups 


2.2 


26 .2 


EC 3.3.-.-: Hydrolases - 
Acting on Ether Bonds 


1.9 


22.2 



CJ 






Protein Family Name 


R- Value 


Value (%) 


EC 2.5.-.-: Transferases - Transferring Alky 1 or 


4.1 


97.2 


Atyl Groups, Other than Methyl Groups 


All lipid-binding proteins 


1.1 


62.2 


Outer membrane 


1.0 


52.6 



ni 






Protein Family Name 


R- Value 


P- Value (%) 


Manganese-binding 


3.5 


96.4 


EC 3. 1 .-.-: Hydrolases - 
Acting on Ester Bonds 


2.9 


93.6 


Iron-binding 


2.1 


25.4 


EC 2.7.-.-: Transferases - 
Transferring Phosphorus- 
Containing Groups 


1.9 


22.2 


All lipid-binding proteins 


1.5 


73.2 


Metal-binding 


1.3 


62.5 



Protein Family Name R- Value P-Value (%) 

All lipid-binding proteins 2.2 26.2 

Metal-binding l.S 80.4 

EC 6.1.-.- : Ligases|- 
Forming Carbon-Oxygen 1.4 
Bonds 

Photosystem I 

Calcium-binding 

TC 3.A.1 ATP-binding 
cassette (ABC) family 

Magnesium-binding 



1.0 
1.0 

1.0 

1.0 



71.3 

52.6 
52.6 

52.6 

52.6 



Fig. 3 Comparative analysis of functional assignment of five different human proteins by SVMProt method 
through BIDD server. A: acetylcholinesterase. B: cytochrome P450. C: gutathione S-transferase. D: protein kinase C. E: alpha-1- 
antichymotrypsin proteins. 
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P450 proteins in forming the pocket are Ile50, Tyr53, 
Argl05, Argl06, Prol07, Glyl09, Glul44, Arg212, 
Phe213, Glu374, and Ile443 (Table 2). Similarly, in 
the case of GST, amino acid residues involved in dif- 
ferent active sites are Tyr6, Trp7, Ile9, ArglO, Arg42, 
Metl04, Asnl08, LeullO, Tyrll5, Aspl61, Lys207 
and Tyr208. The frequently involved amino acid res- 
idues in various active sites of PKC of H. sapiens are 
Leu44, Arg73, Ile74, Ser75, Asp76, Argl22, Glul71, 



Ilel85, Glyl88, Serl89 and Asp202. The last pro- 
tein for analysis is an al-antichymotrypsin which has 
involved residues in different active sites, including 
Asp68, Ser77, Argl23, Asnl25, Vall31, Thrl35, 
Glyl36, Glnl60, Tyrl81, Glnl85, Tyr208, Met282, 
Pro395, and Lys413. The top five different active 
sites predicted through Pocket Finder in five different 
proteins like AchE, P450, GST, PKC and ACT are 
listed in Table 2. 




Fig. 4 The possible binding sites which were obtained by using CASTp server of all five proteins in human. A: 

human acetylcholinesterase (3LII). B: cytochrome P450 (3NXU). C: gutathione S-transferase (4GTU). D: protein kinase C (2XJ1). E: 
a-l-antichymotrypsin (1YXA). The possible binding site indicates in various colors like blue, green, cyan, etc. 



Protein-ligand interaction analysis 

Docking represents the mathematical calculation of 
the most probable spatial orientation of two interact- 
ing molecules, usually protein and small ligand, two 
interacting proteins or DNA and protein. Various pa- 
rameters are calculated to evaluate possibility of such 
protein-ligand interaction. For molecular docking anal- 
ysis, we used new server PatchDock and refinement 
tool FireDock [35] . It has previously been proved that 
all organophosphorus pesticides have toxic effects on 
mammals, birds, fish, reptiles and insects [3] . But there 
is no report about the interaction of these ten organo- 
phosphorus pesticides with target proteins that are in- 
volved in the human metabolic and digestive pathway 
till date. We accessed the tertiary structures from PDB 
with the proposed interaction of five different human 
proteins like AchE (3LII), cytochrome P450 (3NXU), 
GST (4GTU), PKC (2XJ1) and al-antichymotrypsin 



(1YXA) with these organophosphorus pesticides by 
using these tools. They are summarized in Fig. 5 and 
Table 3. The docking solution was visualized using 
the program PyMol V2.0.7 software viewer and dis- 
tance measurements were carried out with the same 
software package. Global energy function (i.e. dock 
score) of five proteins with their highest interaction 
energy with organophosphorus pesticides are calcu- 
lated by the FireDock server to be -45.21 relative units 
(this value is considered to be related to free binding 
energy and higher negative value means higher free 
binding energy and thus higher interaction probabil- 
ity). It is interesting to note that organophosphorus 
pesticide phosmet has shown maximum ligand pro- 
tein interaction (dock score: -45.21) with human P450 
protein in which frequent involvement of H-bonds (4) 
and of both hydrophobic and basic amino acid resi- 
dues including Argl30, Ile443, and Gly444 (Table 4 
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Table 2 The top five involved amino acids residues in active sites in five different human proteins by using 



Pocket Finder server. 



AchE 


P450 


GST 


PKC 


ACT 


SITE:1 


SITE:1 


SITE:1 


SITE:1 


SITE:1 


Gln71, Tyr72 


Ile50, Tyr53, Phe57, Asp76, 


Tyr6, Trp7, Ile9, 


Leu44, Gly45, 


Ser77, Ala80, 


Asp74, Thr 83 


Gln79, Argl05, Argl06, Prol07, 


ArglO, Glyll, 


Gly48, Phe49, 


Leu84, Leul21, 


Trp86, Asn87 


Phel08, Glyl09, Ilell8, Serll9, 


Leul2, Arg42, 


LYS67, Val69, 


Asnl25, Thrl35, 


Pro88, Glyl20, Glyl21, 


Ilel20, Trpl26, Thrl36, Phel37, 


Alal03, Metl04, 


Arg73, Ile74, 


Glyl36, Serl37, 


Glyl22, Tyrl24, Glu202, 


Glul44, Metl45, Thrl87, Serl88, 


Serl07, Asnl08, 


Ser75, Asp76, 


Alal38, Tyrl59, 


Ser203, Trp236, Trp286, 


Lys208, Leu210, Leu211, Arg212, 


LeullO, Alalll, 


Glu89, Lel04, 


Glnl60, Tyrl81, 


Leu289, Ser293,Val294 


Phe213, Ile223, Thr224, Phe241, 


Tyrll5, Leul58, 


Leul20, Argl22, 


Glnl85, Val206, 


Phe295, Phe297 


Arg268, Asp270, Ile301, Phe302, 


Aspl61, Leul65, 


Leul74, Glul71, 


Asn207, Tyr208 


Tyr337, Phe338 


Phe304, Gly306, Glu308, Thr309, 


Hisl66, Lys207, 


Ilel85, Glyl88, 




Tyr341, His447, Gly448, 


Thr310, Arg372, 


Tyr208, Thr209 


Serl89, Asp202 


SITE:2 


Ile451 


Leu373, Glu374, 






Argl23, Leul24, 




Cys442,Gly444 


SITE:2 


SITE:2 


Lysl28, Aspl29, 


SITE:2 


Phe447, Ala448, Leu482, Gln484 


Leu99, Alal03, 


Leul64, Argl66, 


Vall31, Ilel33, 


Pro232, Asn233, Gly234, 




Vall06, Metl33, 


Aspl67, Ilel68, 


Asn269, Phe401, 


Pro235, Trp236, Thr238, 


SITE:2 


Metl34, Phel37, 


Phe201, Tyr207, 


Asp402, Pro409 


Val239, Arg247, Pro290, 


Gly56, Phe57, Cys58, Met59, 


Alal59, Vall62, 


Ser208, Pro209, 




Arg296, Pro368. 


Phe367, Tyr399, Leu475, Leu477, 


Leul63, Hisl66 


Trp212, Tyr218, 


SITE:3 




Leu479, Gly480, Leu482 




Gly220, Ala224, 


Trp215, Lus216, 


SITE:3 




SITE:3 


Ser227 


Val217, Phe219, 


Asn233, Gly234, Pro235, 


SITE:3 


Leu4, Gly5, Tyr6, 




Met244, Leu247, 


Glu313, Val367, Val370, 


Lysl73, Ser312, Ser315, Phe316, 


Argl7, Glu29, 


SITE:3 


Thr249, Leu264, 


His405, Pro410, Gln413, 


Gln484, Pro485, Lys487, Pro488, 


Lys30, Lys31, 


PhelOO, Glyl02, 


Lys265, Leu309 


Trp532, Asn533, Leu536, 


Val489 


Arg201, Leu203 


Glnl50, Vall51, 




Pro537, Leu540 






Ilel73, Leul82, 


SITE:4 




SITE:4 


SITE:4 


Leul84 


Asp68, Cys231, 


SITE:4 


Leul42, Tyr347, Val350, Leu351, 


ArglO, Hisl4, 




Met282, Pro395, 


Arg475, Tyr479, Asn490, 


Arg446, Leu449, Met450, Lys453 


Tyrl60, Aspl61, 


SITE:4 


Lus413, Ala415, 


Glu491, Asp494, Ala497, 




Aspl64, Metl97 


Phel30, Ilel33, 


Asn416 


Pro498, Leu518. 


SITE: 5 




Thrl34, Asp234, 






Asnl59, Aspl74, Alal78, Serl95, 


SITE:5 


Gly238, Asp239 


SITE:5 


SITE:5 


Asnl97, Asnl98, Prol99, Gln200, 


LeullO, Valll3, 




Phe56, Tyr59, 


Glu81, Met85, Aspl31, 


Asp201, Val204, Tyr307. 


Tyrl26, Phel69 


SITE:5 


Leu289, Leu397, 


Vall32, Thr436, Leu437, 






Val96, Ser97, 


Ile398, Met399, 


Trp439, Tyr449, Glu452, 






Ser98,Vall03, 


Ile411, Ala412 


Ile457, Arg463 , Tyr465 






Ilel04, Argl05 





P450: cytochrome P450, GST: gutathione S-transferase, PKC: protein kinase C, ACT: a-l-antichymotrypsin, AchE: acetylcho- 
linesterase. 



and Fig. 6B). The second highest dock score is -39.60 
for GST (Pdbid: 4GTU) with an azinphosmethyl pes- 
ticide as a ligand and the involved hydrophobic amino 
acids, Tyrll5 and Tyr6, are shown in Fig. 6C. AchE 
bound to parathion shows the third highest bind- 
ing energy, which is -37.90 (Table 4 and Fig. 6A), 
and their involved amino acids residues are Tyrl55, 
Tyrl03, Aspl05 and Thrl06, which also belong to the 
active site residues predicted through pocket finder 
and the CASTp server (Table 4 and Fig. 4). Similarly, 
the fourth and fifth highest binding energy is PKC 



bound to azinphosmethyl and ACT bound to quinal- 
phos, which is -36.57 and -30.15, respectively, and the 
number of H-bonding is 4 and 2, respectively. 

In the case of AchE, the binding energy in de- 
scending order with organophosphorus pesticides, 
which exert toxic effects on the human body, varies 
in between -37.90 to -26.53 for parathion, chlorpyri- 
fos, azinphosmethyl, phosmet, quinalphos, fenthion, 
methylparathion, malathion and monocrotophos and 
is shown in Fig. 5. Interaction with the above men- 
tioned pesticides, the frequent involvement of H-bonds 
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Fig. 5 The graphical plot shows the global energy i.e. dock 
score of five different metabolic proteins i.e. acetylcholineste- 
rase (AchE), cytochrome P450 (P450), glutathione S-transferases 
(GST), protein kinase C (PCK) and a-l-antichymotrypsin (ACT) 
with ten different organophosphorous pesticides. 

varies between 1 and 6 and the highest interaction of 
AchE was found to be with methylparathion (6) fol- 
lowed by parathion (5), malathion and monocrotophos 
(4). The amino acids of AchE forming hydrogen bond 
interaction with methylparathion are Trpl82, Argl3, 
and Asnl86, and with parathion are Tyrl55, Tyrl03, 
Aspl05, Thrl06, Serl25 and Ser234, which suggests 
that these pesticides interact with AchE (Fig. 6A) and 
have inhibitory action which occurs as a result of the 
phosphorylation of the serine residue in the active site 
of the enzyme. 

In the case of cytochrome P450n, the interaction 
energy in descending order with several varies in 
between -37.90 to -26.53 is -45.21, -43.97, -43.60, 



-37.23, -36.48, -34.94, -33.09, -31.42 and -25.21, re- 
spectively, for phosmet, azinphosmethyl, quinalphos, 
fenthion, chlorpyrifos, parathion, malathion, methyl- 
parathion, and monocrotophos and the corresponding 
number of H-bonds is 4, 3, 2, 1, 1, 4, 4, 3, and 2, re- 
spectively. The frequently involved amino acids in H- 
bonds are Argl05, Argl06, Arg212, Ala305, Thr310, 
Arg372, Glu374, Ile443, and Gly444 of human P450 
protein (Fig. 6B). P450 isozymes have greater activ- 
ity in direct detoxification of these above mentioned 
-45.21, -43.97, -43.60, -37.23, -36.48, -34.94, -33.09, 
-31.42 and -25.21 respectively; thus, it may be more 
important in assessing risk at low levels of exposure. 

Similarly, in the case of GST, it was reported that 
GST was used for the detection of pesticides like at- 
razine. The nucleophilic attack of GST on atrazine 
releases the H + ion, which can be detected as a pH 
change that directly correlates with the concentra- 
tion of the analyte [11] . The binding energy in descend- 
ing order with various organophosphorus pesticides 
is -39.60, -38.22, and -36.22 for azinphosmethyl, 
phosmet and malathion, respectively, and the lowest 
binding energy is -24.19 for tetrachlorvinphos. For 
interaction of GST with these pesticides, the frequent 
involvement of H-bonds varies between 1 and 4 and 
the highest interaction was found to be with monocro- 
tophos (4) and their dock score is -26.15. Then, mostly 
hydrophobic residues involved in H-bonding are Trp7, 
Tyr6 and Tyrll5. Some other amino acids that are 
involved in the interaction between GST and OPs are 
Arg42, Asnl08 and Lys207. 

The fourth protein, PKC, is also the best target pro- 
tein for organophosphorus pesticides because it can 



Table 3 The global energy (dock score) of five different proteins 
chromeP450 (P450), glutathione S-transferases (GST), protein kinase 
with ten different organophosphorous pesticides 



i.e. acetylcholinesterase (AChE), cyto- 
C (PKC) and a-l-antichymotrypsin (ACT) 



Ligands 


Dock 


No of H bond 


Dock 


No of H bond 


Dock 


No of H bond 


Dock 


No of H bond 


Dock Score 


No of H bond 




Score of 


/Amino acids 


Score 


/Amino acids 


Score 


/Amino acids 


Score of 


/Amino acids 


of 


/Amino acids 




AChE 




of 




of 




PKC 




ACT 










P450 




GST 












Malathion 


-34.30 


4->Tyrl55,Ser234, 


-33.09 


4 -> Prol07, 


-36.22 


2->Asnl08, 


-27.32 


2-> Arg73, 


-23.75 


2->Ser243(2) 






Tyr368(2) 




Argl06(2), Arg372 




TyrllS 




As P 202 






Parathion 


-37.90 


5^Tyrl55, Tyrl03 


-34.94 


4^ Ala305, 


-28.77 


l^Asnl08 


-30.53 


6 -^Serl89, 


-28.04 


2^Gln279 (2) 






(2), Aspl05, Thrl06 




Thr310, Gly444, 








Ser75(2), Arg73 














Ile443 








(3) 






MethylParathion 


-35.22 


6^Trpl32(2), 


-31.42 


3^ Gly444 5 


-28.48 


2 ^TyrllS (2) 


-27.76 


5 -^Arg73, 


-23.35 


4^Serl34(2) 5 






Argl3 (3), Asnl86 




He443, Argl05 








Asp 76(2), 




Glnl32, 


















MetS3(2) 




Tyr210 


Fenthion 


-35.24 


l^Tyrl33 


-37.23 


1^ Arg212 


-32.45 


l->Tyrll5 


-29.89 


l->Ser75 


-24.36 


l->Lys413 


Chlorpyrifos 


-36.23 


2^Tyr337(2) 


-36.48 


1^ Arg212 


-28.46 


2^Trp7(2) 


-32.98 


3^Ser75(2), 


-26.70 


l^Asp68 


















Glyl38 






Monocrotophos 


-26.53 


4^Asp74, Tyr74, 


-25.21 


2^ Ghil25 5 


-26.15 


4-^Trp7, Tjt6 


-21.91 


4^Arg73, 


-16.03 


2^Asp68, 






Tyrl24, Thr75 




Ala289 




(2), Tyrll5 




Phel37, Ser75(2) 




Lys413 


Quinalphos 


36.27 


4->Glu81, Tyr465, 


43.60 


2->Gly431, 


30.66 


2-»Tyr6, 


29.67 


3-»SerlS9, 


30.15 


2->Glu245, 






Arg463(2) 




Argl06 




Lys207 




Asp202 5 GlylSS 




Asp223 


Phosmet 


-36.23 


2^Asn533, Arg417 


-45.21 


4^Argl30 5 Ile443 


-38.22 


l^Tyr6 


-33.23 


5 ^Ser75(3), 


-25.80 


4^Thrl35(2), 










(2), Gly444 








Arg73, Asp 76 




GlnlSS (2) 


Tetrachlorvinphos 


-30.71 


4^Tyel24(2), 


-26.85 


2^ Serl09, 


-24.19 


3^Trp7, 


-29.74 


2 ^Glul71, 


-23.15 


2^Asp68, 






Serl25, Tyr337 




Glu374 




Arg42 (2) 




Asp 178 




Gln283 


Azinpho smethyl 


-36.42 


l^Asp306 


-43.97 


3^Glu374, 


-39.60 


2^Tyrll5, 


-36.57 


4 ^Ser75 (2), 


-25.63 


4^Asnl25, 










Arg372 (2) 




Tyr6 




Asp 76 (2), 




Tyr20S(2), 






















Thrl35 
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Table 4 The docking parameters used for bioinformatics analysis. 



Identified protein's 


Ligands with highest 
docked score 


Global 
energy 


Attractive 
Van der Wall force 


Repulsive 
Van der Wall force 


Atomic contact 
energy (ACE) 


Acetylcholinesterase (Pdbid: 3LII) 


Parathion 


-37.90 


-15.66 


2.83 


-10.50 


Cytochrome protein P450 (Pdbid: 3NXU) 


Phosmet 


-45.21 


-14.83 


2.32 


-15.72 


Gutathione S-transferase (Pdbid: 4GTU) 


Azinphosmethyl 


-39.60 


-15.23 


3.46 


-13.37 


Protein kinase C (Pdbid: 2XJ1) 


Azinphosmethyl 


-36.57 


-13.71 


1.49 


-11.03 


Anti-chymotrypsin (Pdbid: 1YXA) 


Quinalphos 


-30.15 


-14.07 


0.77 


-7.62 



The highest interaction energy of all five proteins like acetylcholinesterase, cytochrome protein P450, gutathione S-transferase, pro- 
tein kinase C and a-l-antichymotrypsin in H. sapiens with the top five different organophosphorous pesticides like parathion, phos- 
met, azinphosmethyl and quinalphos respectively that show the highest interaction with each other. 



also affect signaling pathways by activating PKC [15] . 
However, this activation is probably indirect through 
organophosphorus pesticides-mediated formation of 
reactive oxygen species (ROS) that activate PKC. Or- 
ganophosphorus pesticides that induced ROS forma- 
tion and oxidative stress have also been shown to be 
associated with apoptosis in different tissues. It shows 
that organophosphorus pesticides can inhibit steroid 
androgen receptor, which can cause steroid hormone 
disturbances in the human body. The binding energy 
of PKC with organophosphorus pesticides in descend- 
ing order is -36.57, -33.23, -32.98, -30.53 and -21.91, 
respectively, for azinphosmethyl, phosmet, chlorpy- 
rifos, parathion, and monocrotophos (lowest) (Fig. 5). 
For interaction of PKC with the above mention pes- 
ticides, the frequent involvement of H-bonds varies 
between 1 and 6 and the highest interaction was found 
to be with parathion (6: Serl89, Ser75, and Arg73) 
followed by methylparathion (5: Arg73, Asp76, and 
Met88), phosmet (5: Ser75, Arg73, and Asp76) (Fig. 
6D) that defined their toxicity with PKC. 



In the case of al-antichymotrypsin, which plays an 
important role in the digestive system, the interaction 
energy in descending order with various organophos- 
phorus pesticides is -30.15, -28.04, -26.70, -25.63, 
and -16.03 for quinalphos, parathion, chlorpyrifos and 
monocrotophos, respectively (Fig. 5). The amino ac- 
ids in al-antichymotrypsin frequently involved in the 
formation of hydrogen bond interaction with meth- 
ylparathion are Serl34, Glnl32, and Tyr210, with 
azinphosmethyl, Asnl25, Tyr208, and Thrl35, and 
with phosmet, Thrl35, Glnl85, and Ser234, which in- 
dicates the interaction of organophosphorus pesticides 
with al-antichymotrypsin (Fig. 6E) and organophos- 
phorus pesticides could inhibit the activities of al- 
antichymotrypsin. All the above mentioned five target 
proteins in H. sapiens that are part of the different 
metabolic and digestive pathways are assumed to bind 
with ten organophosphorus pesticides, which inhibit 
the enzymatic activity of these proteins and cause 
toxicity in the human body. 






Fig. 6 Molecular docking through PatchDock and refined by FireDock. A: It represent the acetylcholinesterase (Pdbid: 
3LII) with parathion pesticides and dock score is -37.90. B: The highest binding energy is -45.21 of cytochrome P450 protein (Pd- 
bid: 3NXU) with phosmet. C: Similarly, the glutathione S-tranferase (Pdbid: 4GTU) with ligand as an azinphosmethyl pesticides and 
score is -39.60. D: Binding energy -36.57 of protein kinase C (Pdbid: 2XJ1) with azinphosmethyl. E: Alpha- 1-antichymotrypsin (Pd- 
bid: 1YXA) with quinalphos are -30.15. Structures were visualized using PyMol V2.0.7 software viewer. 
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DISCUSSION 

The crystal structures showed good overall struc- 
tural quality and were validated through PROCHECK, 
ERRAT plot and ProSA program, which revealed 
that PKC indicates a good quality model. Active sites 
and their involved residues are predicted through the 
Pocket finder and the CASTp server. In our dock- 
ing study, the global energy function (i.e. dock score) 
of five proteins with their highest interaction energy 
with organophosphorus pesticides was calculated by 
the FireDock server to be -45.21 relative units. It is 
interesting to note that phosmet pesticides have shown 
a maximum dock score of -45.21 with P450 and the 
involved residues are Argl30, Ile443, and Gly444, 
which also involved in the active sites of the protein 
as shown in Fig. 6B. Thus, our observation suggests 
that in the case of GST bound to azinphosmethyl 
pesticides, the highest dock score is -39.60 and the 
involved hydrophobic amino acids are Tyrll5 and 
Tyr6. AchE bound to parathion shows that the third 
highest binding energy is -37.90 and the involved 
amino acid residues are Tyrl55, Tyrl03, Aspl05, and 
Thrl06, which also belongs to the active sites predict- 
ed through the pocket finder and the CASTp server as 
listed in Table 3 and Fig. 4. Similarly, the fourth and 
fifth highest binding energy was found between PKC 
protein and azinphosmethyl and between ACT and 
quinalphos, which is -36.57 and -30.15, respectively, 
and the number of H-bonding is 4 and 2, respectively. 
We analyzed protein-ligand interaction on five meta- 
bolic proteins interacting with ten organophosphorus 
pesticides and the aim of the study was to identify the 
most toxic organophosphorus pesticides. The dock- 
ing results indicated that phosmet, azinphosmethyl, 
parathion and quinalphos, as the most toxic pesticides, 
inhibit the activity of enzymes and exert adverse ef- 
fects on humans. Further analysis can be carried out 
in the wet lab to determine whether these findings are 
reflective of in vivo conditions. 

In conclusion, the analysis of docking with ten or- 
ganophosphorus pesticides with five target proteins 
(their PDBID are 3LII, 3NXU, 4GTU, 2XJ1 and 
1YXA) in H. sapiens highlighted some important in- 
teractions operating at the molecular level and caused 
toxicity in the human body. The crystal structures 
were shown to possess good overall structural quality 
and were validated by using PROCHECK, ERRAT 
plot and the ProSA program. Active sites and their 
involved residues are predicted through the Pocket 
finder and the CASTp server. The most important as- 
pect of the results is the interaction study in which the 
global energy is maximum at -45.21 for cytochrome 
P450 of H. sapiens with phosmet pesticides, which 



defines the more toxic nature of phosmet that exerts 
harmful effects on humans. AchE bound to parathion 
shows that the third highest binding energy is -37.90. 
The analysis of H-bonding in methylparathion with 
AchE, parathion with PKC and AchE, phosmet with 
PKC and methylparathion with PKC shows that the 
highest number of bond is 6, 6, 5, 5 and 5, respec- 
tively. Therefore, it is concluded that these organo- 
phosphorus pesticides are more toxic and inhibit the 
enzymatic activity by interrupting several metabolic 
and the digestive pathways in H. sapiens. Thus, this 
study will be useful for the prediction of the toxicity 
level of organophosphorus pesticides, herbicides and 
other toxic compounds, which inhibit the activity of 
enzymes. It provides guidance for further screening 
through experimental in vitro and in vivo analyses. 
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